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We test the one-loop chiral perturbation theory formula on unquenched lattice data of pseudoscalar meson decay 
constants. The chiral extrapolation including the effect of the chiral logarithm is attempted and its uncertainty 
is discussed. 



1. Introduction 

In the unquenched QCD simulations using the 
Hybrid Monte Carlo algorithm, the computa- 
tional cost rapidly grows as the chiral limit of 
sea quark is approached. In practical simulations 
with Wilson-type fermions the sea quark mass is 
limited to be heavier than ms/2. To obtain the 
physical results for u and d quarks, therefore, the 
chiral extrapolation is indispensable. 

In the chiral extrapolation the chiral perturba- 
tion theory (ChPT) may be used to decide the 
functional form, as it is an effective theory valid 
for low energy QCD. Once the available lattice 
data are confirmed to be consistent with ChPT, 
the extrapolation to the u and d quark masses us- 
ing the ChPT formula is justified. The practical 
question is, then, whether the lattice results could 
reproduce the sea quark mass dependence, espe- 
cially the chiral logarithm, predicted by ChPT. 

In this talk, we present the sea quark mass de- 
pendence of the pseudoscalar meson decay con- 
stant obtained in unquenched QCD, and compare 
them with the one-loop ChPT prediction. The 
unquenched simulations are done using the stan- 



dard gauge and nonperturbatively 0(a)-improved 
fermion action at (3 = 5.2 with sea quark masses 
corresponding to the pion mass in the range 550- 
1000 MeV. Further details of the simulation are 
discussed in a separate talk 

We also discuss the uncertainty associated with 
the chiral extrapolation taking light-light and 
heavy-light decay constants as examples. 

2. Test of the ChPT formula 

In full QCD the ChPT predicts a specific func- 
tional dependence of physical quantities on the 
quark mass, i.e. the chiral logarithm, at the one- 
loop order. For Nf flavors of degenerate quarks 
with a mass ms, the pseudoscalar meson decay 
constant /ss is given by |^ 
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with yss — '^Bqius / {^11 fY . While the low energy 
constants are unknown parameters, the chiral 
log term yss In 2/ss appears with a definite coeffi- 
cient depending only on the number of flavors. 

Figure [l] shows the lattice results together with 
fitting curves. If we leave the coefficient of the 
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Figure 1. Fit of pseudoscalar meson decay con- 
stant with the one-loop ChPT formula. 
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Figure 2. Ratio test of the partially quenched 
ChPT formula. 



chiral log term as a free parameter, the fit result 
is consistent with zero (solid line), while the fit 
with the fixed coefRcient Nf/2 gives a bad x^/dof 
(dashed line). A similar observation is obtained 
for the PCAC relation M^g/2ms §. 

Partially quenched ChPT may be used 

to explicitly test the presence of the chiral log- 
arithm. For non-degenerate pions composed of 
quarks with mass ms and my {V denotes a va- 
lence quark) , the low energy constants cancel out 
in the double ratios 
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with t = \n{yvv /yss) + 1 - (yvv/yss), and only 
the chiral log terms remain. A parameter yss/Nf 
obtained as a coefficient of t in (H) and in ( |) are 
plotted as a function of Mgg in Figure |. We 
find that the results are much smaller than the 
prediction of the partially quenched ChPT shown 
by a steep dashed line. 

3. Uncertainty in the chiral extrapolation 

Since the lattice data do not support the pres- 
ence of the chiral logarithm for the sea quark 
masses used in the simulation, the chiral extrap- 
olation using the one-loop ChPT formula is not 



fully justified. Instead, we consider several possi- 
ble functional forms to approach the chiral limit 
and discuss their associated uncertainty. 

Our observation suggests that the mass region 
where the chiral logarithm becomes important is 
around or below 500 MeV, and the ChPT ceases 
to converge above that scale. Then, a possible 
way to extrapolate the data including the effect 
of the chiral logarithm is to use a polynomial fit 
(quadratic fit, for example) above some energy 
scale /i, and then switch to the one-loop ChPT 
formula below ^. An example is shown in Figure ^ 
for ^l = 300 and 500 MeV. The limit of /i = MeV 
corresponds to the usual polynomial fit. Since the 
scale n is unknown, the variation of several fit 
curves, about ± 10% in the chiral limit, should 
be taken as systematic uncertainty. 

Another possible functional form suggested by 
the Adelaide-MIT group § is the one-loop ChPT 
with a hard momentum cutoff /i. It amounts 
to replace the chiral log term ln(mj/^^) by 
TOjln(m^/(m^ + /i^)). Changing the unknown 
"cutoff" scale fj, from to 1 GeV, we obtain the 
similar size of uncertainty in the chiral limit. 

The SU(3) breaking ratio in the decay con- 
stant /aV/tt niay be obtained with a fit to par- 
tially quenched lattice data. Our result from 
the above fit functions changes from 1.167(3) 
(quadratic fit, corresponding to /i = MeV) or 
1.190(3) (Adelaide-MIT fit, ^JL = 500 MeV) to 
1.276(7) (chiral log plus quadratic, correspond- 
ing to /.t = CX3 MeV). Although it is clear that our 
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Figure 3. Uncertainty in the chiral extrapola- 
tion. 

two-flavor QCD result is significantly higher than 
the quenched result f .081(5)(17) § and closer to 
the physical value 1.22, the uncertainty is still 
sizable. 

For the heavy-light decay constant the predic- 
tion of ChPT is available in the heavy quark limit 
for quenched, partially quenched and full QCD 
The chiral logarithm appears with a def- 
inite coefficient but including an additional cou- 
pling constant g describing the B* Bit interaction. 
The chiral extrapolation of fs and fs^ including 
the chiral logarithm is shown in Figure ^ with 
two representative values oi g. As in the pion de- 
cay constant, the uncertainty in the chiral limit 
is enhanced by the chiral logarithm. 

The ratio Jb^ / fs is needed in the extraction 
of the CKM matrix element |Vtrf/Vts|. Since the 
bulk of systematic errors cancels in the ratio, one 
may expect better accuracy than the determina- 
tion of \Vtd\ solely from AM^. Our preliminary 
result varies from 1.24 (quadratic fit) to 1.38 (chi- 
ral log, g — 0.59). It suggests that the ratio and 
its error can be significanly larger than the pre- 
vious world average 1.16(4). Similar discussion, 
but using quenched data, has been made in [l0| . 

4. Conclusions 

The chiral logarithm expected from ChPT is 
not observed in the unquenched lattice data with 
rrips greater than about 500 MeV, which suggests 
that ChPT may only be applied in smaller mass 
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Figure 4. Chiral extrapolation of heavy-light 
decay constants. 

regions. The estimate using model functions for 
the chiral extrapolation leads to the uncertainty 
as large as ±10% for the decay constants. 
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